Epidemic Spreading in Weighted Networks: An Edge-Based Mean-Field Solution 
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Weight distribution largely impacts the epidemic spreading taking place on top of networks. This 
paper studies a susceptible-infected-susceptible model on regular random networks with different 
kinds of weight distributions. Simulation results show that the more homogeneous weight distri- 
bution leads to higher epidemic prevalence, which, unfortunately, could not be captured by the 
traditional mean- field approximation. This paper gives an edge-based mean- field solution for gen- 
eral weight distribution, which can quantitatively reproduce the simulation results. This method 
could find its applications in characterizing the non-equilibrium steady states of dynamical processes 
on weighted networks. 

PACS numbers: 89.75.Hc, 05.45.Xt, 05.70.Ln, 87.23.Ge 



I. INTRODUCTION 



Many social, biological and communication systems 
can be described as networks with nodes denoting indi- 
viduals and edges representing interactions among them. 
Edges in many networks do not only indicate the ex- 
istence of interactions, but also have associated with 
weights to record the interacting strengths. Examples 
are numerous, including the scientific collaboration net- 
works where the edge weight between two scientists de- 
pends on the number of co-authorized papers and the 
number of authors in each of these papers the air- 
port networks with edge weight representing the number 
of flights Q or the number of available seats @ between 
two airports, the food webs where edge weight can be 
defined by interaction frequency, by carbon flow or by 
interaction strength [4j, the cell-phone communication 
networks with weight being assigned as the total number 
of calls or the total duration of calls [E| , the world trade 
web where the edge weight accounts for the exports over 
GDP of the exporting country [6|, and so on. 

Weight distribution has remarkable effects on the prop- 
erties of dynamical processes taking place on top of net- 
works, such as synchronization [7H9|, transportation and 
routing [Iol4l2L evolutionary game [l3[ [liL cascading 
failu re Il5l Il6j , diffusion [17], [H| , control [19[ , condensa- 
tion [20], and so forth. Epidemic spreading in weighted 
networks has also been considered in the literature [231 ] . 
Yan et al. [24[ studied a susceptible-infected (SI) model 
on weighted scale-free networks and found that the het- 
erogeneous weight distribution will slow down the spread- 
ing. Therefore, given the topology and mean infectiv- 
ity, the epidemic spreads fastest in unweighted networks. 
Chu et al. [25| investigated the same dynamical pro- 
cess upon weighted scale-free networks with community 
structure, and showed that the weights of edges among 
different communities have higher im pac ts than those of 
edges inside communities. Chu et al. [26j further studied 
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the susceptible-infected-recovered (SIR) model and found 
that the weight distribution has considerable impacts on 
both epidemic threshold and prevalence. Karsai et al. 
p7j considered the contact process in weighted scale-free 
networks, in which the weight of an edge connecting two 
higher-degree nodes is relatively small. Yang et al. [Hf 
further proved that in the contact process, the epidemic 
prevalence can be maximized by setting the egde weight 
inversely proportional to the degree of receiving node. 

Mean-field approximation is a powerful tool to ana- 
lyze dynamical processes on networks. In particular, the 
so-called heterogeneous mean-field approximation has 
achieved great success in predicting threshold of epidemic 
spreading in scale- free networks jlH, [30[. However, this 
method may fail to capture the behavior of some dy- 
namical processes in the thermodynamic limit (3l| and 
can not account for the fluctuations caused by quenched 
structure |32j |. In this paper, we study epidemic spread- 
ing in homogeneous networks with heterogeneous weight 
distribution, and develop an edge-based mean-field ap- 
proximation. Simulation results indicate that the more 
homogeneous the weight distribution is, the higher the 
epidemic prevalence is. This phenomenon can be quan- 
titatively reproduced by the edge-based mean-field solu- 
tion. In contrast, the node-based mean- field solution is 
insensitive to the change of weight distribution, and it 
tends to overestimate the epidemic prevalence. 

This paper is organized as follows. In the next section, 
the model will be introduced, including how to construct 
networks, how to generate weight distributions, as well 
as the details of spreading dynamics. In Section 3, the 
edge-based mean-field solution will be presented, while 
for comparison, the node-based mean-field solution will 
be introduced in the Appendix. Section 4 compares the 
simulation results and analytical solutions and give an ex- 
planation about why the edge-based mean-field solution 
is more accurate than the node-based one. Eventually, 
we draw the conclusion, emphasize the advantages of the 
present method and discuss its potential relevance and 
applications in the last section. 
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II. MODEL 

To emphasize the effects caused by the weight distri- 
bution, we concentrate on the regular random networks 
where edges are added in a completely random matter 
and the degree of every node is fixed as k. In this way, 
the impacts of local clustering and heterogeneity in net- 
work topology can be eliminated. 

We investigate two sets of weight distributions. In the 
former case, the weight of each edge is the average of n 
elements independently drawn from the uniform distribu- 
tion U(0, 2). When n = 1, the weight distribution p(w) is 
exactly the same as U(0, 2), and p(w) will become more 
homogeneous with the increasing of n, following the form 



J>K1) = U(0,2), 

. f 2 ,nw — z 

p(w;n) = / zp( ;n — l)az,n = 2,3, 

Jo n-1 

a 2 



(1) 



lim p(w; n) ~ N(/i, — ), 



where N stands for the normal distribution, and \i = 1 
and a = yj x 2 dx « 0.57735 are the mean value and 
standard deviation of the uniform distribution U(0,2), 
respectively. Figure 1(a) illustrates example distributions 
for n = 1, 4, 7 and 10. 

In the latter case, the edge weight distribution follows 
a power law 



p(w) 



(2) 



The power-law distributions are generated according to 
the method in Refs. [U, [22| with weights between 
^min = 1 and i^max = 10 4 . Figure 1(b) illustrates ex- 
ample distributions with 7 = 2, 3 and 4. 

The epidemic spreading is described by a weighted 
susceptible-infected-susceptible (SIS) model. Each node 
is in one of two discrete states, either susceptible (S) or 
infected (I). At each time step, each infected node will 
select one of its neighbors to contact with. The selection 
probability is proportional to the edge weight. If the se- 
lected neighbor has been infected already, nothing hap- 
pens, while if it is susceptible, with probability A, it will 
be infected and then turn to be I-node in the next time 
step. Meanwhile, each infected node will become suscep- 
tible in the next time step with probability f3. We assume 
that each infected node can contact only one neighbor per 
time step since in real social networks even an individual 
may have many friends, her could only interact with a 
very limited number of them whatever by physical con- 
tacts or by information communication [33|, [34| . 



III. ANALYSIS 

The analytical results can be obtained by the Markov 
chain equations about the evolution of edge states. For 





FIG. 1: (Color online) Illustration of weight distributions: (a) 
more homogeneous distributions, as described in Eq. 1 and (b) 
more heterogeneous distributions, as described in Eq. 2. Each 
distribution is obtained by sampling 10 6 weights. 



every edge in the network, there are three possible states 
depending on its two ends, say SS, SI and II. With an 
edge-based mean- field hypothesis, we assume that edges 
with the same weight distribute uniformly in the network. 
Given a specific weight w, we denote ess( w ), ^si( w ) an d 
eu(w) as the densities of edges in corresponding states, 
which satisfy the normalized condition 



ess(w) + e S i(w) + e n (w) 1. 



(3) 



In the steady state, not only the density of the 
infected node, but also the densities of edges in 
different states keep stable. Denote w) = 

[ess(t, w), esi(t, w), esi(t, u))] T the densities of edges in 
different states at time step t. A time-dependent trans- 
ferring matrix, T(t,w), governs how 1^(t,w) changes to 
, w) during the time step £, so the Markov chain 
equations can be written as, 

~^(t + 1, w) = T(t, w)^(t, w). (4) 
An infected node, no matter it associates with which 



kinds of edges, will change to a susceptible one with prob- 
ability j3. Clearly, the more infected nodes a susceptible 
node connects, the higher probability it will be infected. 
However, we could not precisely know each susceptible 
node's local surrounding, so we only consider two dif- 
ferent prior conditions: a susceptible node is one end 
of an Si-edge or it is one end of an SS-edge. We de- 
note (p(t) the probability that a susceptible node con- 
necting to an SS-edge with weight w will become an in- 
fected node during the time step £, then for a suscep- 
tible node connecting to an Si-edge with weight w, the 
probability it will be infected during the time step t is 



i - [i - <p{t)\ 

weight. The item 



{k-l)(w)+w 



where (w) is the mean 



1 



-A 



is the probability 



(k-l)(w)+w ' 

that the monitored susceptible node will not be infected 
by the known infected neighbor with Tl — ,w v , — denot- 

J to (k— l)(w)+w 

ing the probability this infected neighbor will choose the 
monitored node to contact with, meanwhile the item 
[1 — (p(t)] is the probability that the monitored node will 
not be infected by the other k — 1 neighbors. Notice that, 
the item ( k _ 1 y^ w ^+ w is an approximation only valid for 
homogeneous weight distribution, which is good enough 
for the case shown in figure 1(a) yet not quantitatively 
suitable for the case shown in figure 1(b), which will be 
discussed later. 

Denoting A, B and C the probabilities that during the 
time step t, an infected node changes to a susceptible one, 
an susceptible node as one end of an SS-edge weighted w 
changes to an infected one and an susceptible node as one 
end of an Si-edge weighted w changes to an infected one, 
respectively, then B and C depends on t and w. Without 
lose of generality, we always adopt a normalized weight 
distribution with (w) = 1, and the changing probabilities 
can be summarized in the follows. 



A = /3, 
B = <p{t), 

C = 1 - [1 - <p(t)] ( 1 



k - 1 



■ w 



(5) 



We neglect the indices t and w of B and C for simplic- 
ity. Accordingly, the transferring matrix T(£, w) for the 
changes of densities of edges in different states is 



(1-B) 2 
2(1 — B)B 
B 2 



A(l - C) A 2 
A-C + 2AC 2A(1 — A) 
(1-A)C (1-A) 2 



(6) 

As a transferring matrix, T is column normalized. 

Under the above framework, the estimation of cp(t) be- 
comes the most significant part. Each susceptible node 
is associated with k edges, which can only be in SI or SS 
states. For example, as shown in figure 2, the suscepti- 
ble node in the center is surrounded by k = 10 nodes, 
and the density of its associated Si-edges is 0.5. Denote 
by psi the average density of Si-edges associated with a 
randomly selected susceptible node, if the edges in differ- 




FIG. 2: (Color online) Illustrations of the densities of Si-edges 
and active Si-edges. The central node is a susceptible node 
with 10 neighbors, among which, five are susceptible nodes 
and the other five are infected nodes. Therefore, the density 
of Si-edges for the central node is 0.5, and psi is the average 
density of Si-edges over all susceptible nodes. Among these 
five infected neighbors, only two of them choose the central 
node to contact (indicated by red color and thick lines), there- 
fore the density of active Si-edges is 0.2. Analogously, p is the 
average density of active Si-edges over all susceptible nodes. 



ent states distribute uniformly in the network, the num- 
ber of Si-edges and SS-edges are Msi = Nskpsi and 
Mss — \Nsk{\ — psi), where Ns denotes the number 
of S-nodes and the factor 1/2 accounting for the dupli- 
cate count of SS-edges. Therefore, the density psi can 
be written as 



Psi 



Msi 



2M SS + Msi 



(7) 



This simple derivation could find its application in esti- 
mating (f(t) by further considering (i) edges with differ- 
ent weights should be treated separately to ensure the 
accuracy, and (ii) an Si-edge is active only if the cor- 
responding /-node chooses this edge to make a contact. 
For example, in figure 2, the density of Si-edges subject 
to the central node is only 0.2. Accordingly, the average 
density of active Si-edges associated with a randomly se- 
lected susceptible node is 



P(t) 



J e S i(t,w)p(w) 



k—l+w 



dw 



f[e S i(t, w) + 2e S s(t, w)]p(w)dw ' 



(8) 



where is the approximate probability an Si-edge 

with weight w gets active with homogeneous weight dis- 
tribution. If there are I active Si-edges connecting to 
a susceptible node, it will be infected with probability 
1 — (1 — X) 1 . With the density of active Si-edges being 
p(t) and the degree of each node being fc, the infected 




300 



k-i+w •> ^ s no ^ accura ^ e enough. In this case, we 
replace it by 



(j) + w 



q{4>)d<i), 



(11) 



where </> denotes the sum of k — 1 randomly selected 
weights from the distribution p(w) and q(<p) is the distri- 
bution of cj). Generally speaking, given p(w), q((j>) can be 
easily determined numerically. Therefore, via replacing 



the term 



in Eq. 5 and Eq. 8 by Eq. 11, we can 



k— l+u> 

more accurately characterize the systems with heteroge- 
neous weight distributions. 




FIG. 3: (Color online) Density of infected nodes versus time 
in regular random networks with N = 10 3 and k = 5. Ini- 
tially, 10 nodes are randomly selected to be infected, and the 
infecting rate and recovering rate are A = 0.2 and f3 = 0.1, 
respectively. Every data point is averaged over 10 3 indepen- 
dent realizations. Subgraphs (a) and (b) imply the spreading 
processes in the networks with different kinds of weight dis- 
tributions. Three curves in subgraph (a) correspond to Eq. 1 
with parameter n = 1, 5 and 10, while those in subgraph (b) 
to Eq. 2 with parameter 7 = 2.0, 3.0 and 4.0. 

rate of a susceptible node in an SS-edge is 
ip{t) = ^C l k _ lP l (t)[l -pit)}*- 1 - 1 [1 - (1 - A)'] . (9) 

By iteratively multiplied by T(t,w), w) will 

rapidly converge to a steady vector E(w) = 
[ess(w),esi(w),eu(w)] T , and the density of infected 
nodes in the non-equilibrium steady state can thus be 
obtained as 

i = \ J i 2e id w ) + e S i(w)]p(w)dw. (10) 

If the weight distribution p(w) is highly heterogeneous 
like those shown in Fig. 1(b), the estimation of the prob- 
ability an infected node will choose an edge with weight 



IV. RESULTS 

We have checked that the steady state is insensitive 
to the initial condition and the network size, and thus 
initially we always randomly select 10 nodes as infected 
and all the shown results are implemented on a network 
of 10 3 nodes and averaged over 10 3 independent realiza- 
tions. Figure 3 displays the fundamental graph of an 
SIS model, say the density of infected nodes, i(t), as a 
function of time. Note that, p(w) will become more ho- 
mogeneous as the increasing of n and 7, and at the same 
time, the density of infected nodes in the steady state 
increases. In other words, the more homogeneous weight 
distribution leads to higher epidemic prevalence. 

To further explore how the weight distribution affects 
the density of infected nodes in the steady state, z, we 
show in figure 4 how i changes with n and 7. It is consis- 
tent with figure 2 that the larger n or 7 is, the more in- 
fected nodes exist in the stable state. In order to empha- 
size the advantages of the edge-based mean-field solution, 
not only do we compare it with simulation results, but 
with the classical node-based mean-field solution as well, 
which is introduced in Appendix A. Very clearly, the 
edge-based mean-field solution can quantitatively repro- 
duce the simulation results, while the node-based mean- 
field solution completely fail to capture the impacts of 
weight distribution. 

Figure 5 provides an evidence that could, to some ex- 
tent, explain the advantage of edge-based mean- field ap- 
proximation. Even for the case n = 1 corresponding to 
the most homogeneous weight distribution discussed in 
this article, the densities of edge states change with the 
weight w. This will become even more remarkable when 
the weight distribution gets broader. The edge-based 
mean-field solution (see figure 5b) can quantitatively cap- 
ture the weight-dependent state densities (see figure 5a). 
This figure also clearly points out the causation of the 
failure of node-based mean-field approximation, that is, 
the ignorance of the different density of edge states at 
different weights. 
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FIG. 4: (Color online) How the density of infected nodes in 
the stable state changes with the weight distribution with 
N = 10 3 , A = 0.2 and f3 = 0.1 fixed. The red and blue 
colors represent the cases k = 5 and k = 10, respectively. 
Simulation results (triangles and circles) are the average over 
10 3 independent realizations, with the error bars denoting the 
standard deviations. The solid and dashed lines correspond to 
edge-based and node-based mean- field solutions, respectively. 
Similar to figure 2, subgraphs (a) and (b) are for different 
kinds of weight distributions. 



V. DISCUSSION 

To solve the SIS model on weighted regular random 
networks, this paper proposed an edge-based mean- field 
method, which outperforms the traditional node-based 
method and can quantitatively reproduce the simulation 
results. Both analytical solution and simulation results 
show that the more homogeneous weight distribution 
leads to higher epidemic prevalence. 

The main shortcoming of the mean- field approximation 
in network analysis is its disability in capturing the fluc- 
tuations caused by quenched structure. Compared with 
the traditional node-based mean- field method [29|, l30| , 
more accurate solutions could possibly be obtained by 
extending the analyzed units from nodes to edges, form 
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FIG. 5: (Color online) How the densities of edges in different 
states (II, SI and SS) change with the weight w in the steady 
state. The weight distribution is set as U(0, 2) (i.e., n = 1) 
and other parameters are N = 10 3 , k = 10, A = 0.2 and 
P — 0.1. Subgraphs (a) and (b) are the simulation result 
(averaged over 10 3 independent realizations) and analytical 
solution, respectively. 



edges to 3-order subgraphs, and so on. The mean-field 
approximation for correlated networks [35|, to some ex- 
tent, can be considered as an elegant example in this 
direction since the degrees of two ends of an edge are si- 
multaneously considered. The present work also goes in 
such direction, but it directly considers the heterogeneity 
of edge weights instead of that of node degree. 

The present method could definitely find applications 
in analyzing dynamical processes on weighted networks, 
such as cascading, voting model, evolutionary game, 
transportation, etc. However, in reality, both node de- 
grees and edge weights can be highly heterogeneous, and 
neither of the traditional node-based mean-field approx- 
imation nor the present edge-based mean-field approxi- 
mation is able to give very accurate predictions. How to 
integrate these two methods is still an open issue. 
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Appendix A: Node-Based Mean-Field Solution 

Denoting N$(t) the number of susceptible individuals 
at time step t and psiif) the average density of Si-edges 
associated with a randomly chosen susceptible node, the 
number of Si-edges and SS-edges are 



Mssit) = \N s (t)k(l - psi(t)), 
M S i(t) = N s (t)kpsi(t). 

Therefore, the density psi(t) can be written as 

M S i(t) 



(Al) 



Psi(t) = 



M SI (t) + 2M ss (t)' 



(A2) 



Denoting i(t) the density of infected individuals, accord- 
ing to the mean-field approximation, psi(t) = ~ 
i(t) when N ^> 1. Similar to the edge-based analysis, 
the average density of active Si-edges associated with a 
randomly selected susceptible node is 



p(t) = p S i(t) 



J pH J 

i(t) j p(w) J ■ 



w 



<p + w 
w 



q((j))d(j)dw, 



(A3) 



q((j))d(j)dw. 



(j) + W 

Therefore, the infected rate of a susceptible node is 

<p(t) = tf-lC l k _ lP (t) l (l - p®)"- 1 -^ - (1 - A)']. (A4) 

Knowing the infecting rate and the recovering rate, /3, 
the density of infected individuals evolves as 

i(t + 1) = i(*)(l - /?) + [1 - i(t)] ¥>(*). (A5) 

It will rapidly converge to a stable value i. The results 
are shown as dashed lines in figure 3. 
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